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Abstract 

Using an efficient computational approach, we have reconstructed the structure of 
the dust cloud in the Solar system between 0.5 and 100 AU produced by the Kuiper belt 
objects. Our simulations offer a 3-D physical model of the 'kuiperoidal' dust cloud based 
on the distribution of 280 dust particle trajectories produced by 100 known Kuiper belt 
objects ; the resulting 3-D grid consists of 1.9 x 10^ cells containing 1.2 x 10^^ particle 
positions. The following processes that influence the dust particle dynamics are taken into 
account: 1) gravitational scattering on the eight planets (neglecting Pluto); 2) planetary 
resonances; 3) radiation pressure; and 4) the Poynting-Robertson (P-R) and solar wind 
drags. We find the dust distribution highly non-uniform: there is a minimum in the 
kuiperoidal dust between Mars and Jupiter, after which both the column and number 
densities of kuiperoidal dust sharply increase with heliocentric distance between 5 and 10 
AU, and then form a plateau between 10 and 50 AU. Between 25 and 45 AU, there is an 
appreciable concentration of kuiperoidal dust in the form of a broad belt of mostly resonant 
particles associated with Neptune. In fact, each giant planet possesses its own circumsolar 
dust belt consisting of both resonant and gravitationally scattered particles. As with the 
cometary belts simulated in our related papers (Ozernoy, Gorkavyi, & Taidakova 2000a,b), 
we reveal a rich and sophisticated resonant structure of the dust belts containing families 
of resonant peaks and gaps. An important result is that both the column and number dust 
density are more or less flat between 10 and 50 AU, which might explain the surprising 
data obtained by Pioneers 10 & 11 and Voyager that the dust number density remains 
approximately distance-independent in this region. The simulated kuiperoidal dust, in 
addition to asteroidal and cometary dust, might represent a third possible source of the 
zodiacal light in the Solar system. 

Keywords: Interplanetary dust, the Edgeworth-Kuiper belt 
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1. Introduction 

It is well known tliat tlie Edgewortli-Kuiper belt objects (called 'kuiperoids' hereinafter) 
can replenish the cometary populations throughout the Solar system (e.g., Levison & 
Duncan 1997; Ozernoy, Gorkavyi, & Taidakova 2000a,b = OGT 2000a,b; and refs. therein). 
Recently, it has been recognized that kuiperoids might be also one of the major sources 
of dust in the Solar system (e.g., Backman et al. 1995, Liou et al. 1996). This dust can 
be produced due to evaporation of the volatile material from the surface of kuiperoids by 
processes, such as the solar radiation and wind, mutual collisions of kuiperoids, micrometeor 
bombardment, etc. In this paper, we perform extensive numerical simulations to examine 
the distributions (both in orbital parameters and in space) of kuiperoidal dust particles and 
thereby to analyse the structure of the kuiperoidal dust cloud. 

The dynamics of this dust is determined by three major effects: (i) the Poynting- 
Robertson (P-R) drag (including radiation pressure and solar wind drag), (ii) gravitational 
scattering on the planets, and (iii) resonances with the planets. Extensive work on dust 
particle evolution governed by these effects has been done by a number of investigators 
(WeidenschiUing & Jackson 1993; Hamilton 1994; Roques et al. 1994; Liou & Zook 1997; 
Gorkavyi, Ozernoy, Mather, & Taidakova (= GOMT) 1997a,b and 1998a,b; Kortenkamp & 
Dermott 1998). The evolution of kuiperoidal dust was analyzed by Liou, Zook, & Dermott 
(1997), Gorkavyi, Ozernoy, & Taidakova (1998), and Liou & Zook (1999). The present 
paper makes a next step by employing a fast and efficient method to compute a stationary 
distribution function of dust particles in the phase space, which provides much better 
statistics to derive a 3-D model of the interplanetary dust cloud and to reveal its rich 
resonant structure. 

In Sect. 2, we discuss the sources of dust in the outer Solar system. The dynamical 
evolution of dust particles is reviewed in Sect. 3. Sect. 4 describes our numerical method 
that enables us to compute the 3-D distribution of dust in the Solar system. We employ 
an implicit second-order integrator for dissipative systems (Taidakova & Gorkavyi 1999) 
outlined in Appendix. Sect. 5 contains the results of these computations, which reveal 
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the global dust distribution as well as interesting details of its resonant structure. Our 
conclusions are presented in Sect. 6. 

2. The Sources of Dust Particles 

There is mounting evidence that the sources of the interplanetary dust particles 
(IDPs) cannot be entirely reduced simply to those comets which produce the observed 
dust tails and/or to asteroids which are thought to be responsible for the observed 'dust 
bands' in the IDP emission. Many facts forces us to suspect that additional sources of the 
interplanetary dust must exist. Among others, two groups of facts are worth mentioning: 
(i) According to Pioneers 10 & 11 and Voyagers 1 & 2 data, the dust number density is 
approximately distance-independent in the outer Solar system between 10 and 40-50 AU 
(Humes 1980, Divine 1993, Gurnett et al. 1997), while both the asteroidal and cometary 
dust number densities are known to sharply decrease with heliocentric distance (GOMT 
1997b); (ii) Chemical analyses and other space-based data indicate that some IDPs spent a 
much larger time in space than the typical asteroidal and cometary particles (Flynn 1996). 
Thus, there is strong evidence in favor of other sources of dust in the Solar system, along 
with the known comets and asteroids. 

This third component of the IDP cloud might be the 'kuiperoidal' dust (Backman et 
al. 1995). The Kuiper belt can influence the formation of the IDP cloud in two ways: 
1) as a source of small-size particles slowly drifting toward the Sun under the combined 
action of the PR-drag, gravitational scattering, and resonances; and 2) as a source of 
trans- Jovian comets. It is commonly agreed that the Jupiter-family comets are produced 
by transporting the comets from the Kuiper belt via gravitational scattering on the four 
giant planets so that each planet scatters the comets both toward and away from the Sun 
(Levison & Duncan 1997). Our simulations (OGT 2000a,b) indicate that, between Jupiter 
and Neptune, there is a large population of minor bodies forming four cometary-asteroidal 
belts near the orbits of the giant planets. 

According to our simulations, the minor body families of Saturn, Uranus, and Neptune 
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should contain progressively larger numbers of comets than one sees near Jupiter. Even 
despite a many-fold decrease of the solar heat intensity at such large distances, those 
numerous comets may produce dust in amounts comparable to that from a few active 
J-comets. Complementary mechanisms of dust release from kuiperoids and Centaurs 
between Jupiter and Neptune can include impacts of large grains and the solar wind. 
Without discussing here the dust production by the above mechanisms, which is out of 
our paper's scope, we simply refer to observational data that indicate, for a number of 
kuiperoids and Centaurs, a steady cometary activity lasting for years (e.g. Brown & Luu 
1998 and refs. therein). 



The kuiperoidal dust experiences the same dynamical effects as the cometary and 
asteroidal dust, with the only difference that, due to a slower PR-drift and a stronger 
influence of the giant planets, the role of gravitational scattering and resonance captures 
must be more important for it. 

Just after the birth of a dust particle, the solar pressure forces it to change its orbit to 
a more distant and eccentric one. This change of the orbital parameters (semimajor axis 
ttd and eccentricity ea) of dust particles, which start their journey from the apocenter or 
pericenter of kuiperoidal orbits, is described by 



Here qk and ck are the orbital parameters of kuiperoids, signs and "— " correspond to 
the start from apocenter and pericenter, accordingly, and P LQ/{MQr) is the ratio of the 
solar light pressure and the gravitational force applied to the dust grain of radius r (in /xm). 
In what follows, we consider two /9-values: (3 — 0.285 and (3 — 0.057, which correspond, for 
grains of density p— {1 — 2) g/cm^, to r = (1 — 2) //m and r = (5 — 10) //m, respectively. 

The basic mechanisms of the dynamical evolution of dust are as follows: 



3. Dynamical evolution of dust particles 



l-f3 



(1) 
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1. Drift of particles toward the Sun under the Poynting- Robertson drag. Such drift 
occurs more or less freely except for short periods when particles experience a strong 
gravitational scattering on the planets or when they pass through the outer or inner 
resonances without being captured into them. A particle governed by the P-R drag drifts 
in the (a, e, i)-space, with a decreasing a and e, along the line 

a(l — e^)e~^^^ = const, i = const. (3) 

2. Gravitational scattering by the planets, which builds up a population of particles 
with large eccentricities, including grains that are ejected from the Solar system. A 
gravitationally scattered particle experiences a rather chaotic 'jump' in the (a, e, i)-space, 
but conserves its Tisserand parameter T during the jump: 

Although the P-R drag changes that parameter, it happens on a much longer timescale. 
As shown in OGT (2000b), gravitational scattering results in the motion of most of the 
particles within the so called 'crossing zone' (i.e. the zone of strong gravitational scattering) 
defined, in the (a, e)-plane of orbital coordinates, by: 

a(l — e) < api if a > api, 

a(l + e) > ttpi if a < Gpi. (5) 

Here a and e are the semimajor axis and eccentricity of the test particle, respectively; and 
ttpi is the semi-major axis of the planet. As a result of gravitational scatterings, the particle 
pericenters are close enough to the planet's orbit Gpi (OGT 2000b): 

a(l — e) ~ Gpi if a > Gpi. (6) 

3. Resonant capture of dust into outer resonances with the planets. A particle captured 
into a resonance is positioned on the line 

a^il-f3y/'[^Y\pU (7) 
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where usually n > m (i.e. this is an outer resonance). The eccentricity of the resonant 
particle increases with time (Weidenschilling and Jackson 1993), while its i oscillates and 
gradually decreases (Liou & Zook 1997). 

Athough the above dynamical trajectories governed by just one dominating dynamical 
factor turn out to be only approximate as soon as the two other factors are taken into 
account, they are very helpful as analytical approximations in the appropriate limiting 
cases. Moreover, in accordance with the dominating dynamical factor, we expect to get 
three major components of dust populations: i) 'freely' drifting particles, (ii) gravitationally 
scattered particles, and (iii) particles captured into resonances. This classification is helpful 
while interpreting the results of the present numerical simulations. In Sect. 5, wc consider 
the steady-state distribution of these three dust components, but before that we describe 
our computational approach. 

4. Computational Method: Simulation of a Quasi-statiouciry Distribution of 

Dust Particles in the Solar System 

We calculate the orbital elements (a, e, i) and then the spatial positions of massless 
particles starting from a particular kuiperoid as the source of dust and drifting toward the 
Sun under the P-R drag. On its way to the Sun, each particle undergoes the gravitational 
influence of the planets. To save computational effort, we assume that the Sun is flxed at 
the origin and the 8 planets (excluding Pluto) are on circular orbits with zero inclinations 
(this approximation will be abandoned in our further work). 

The distribution of dust in the Solar system averaged over a time scale of planetary 
orbital motions is described as stationary. To simulate this stationary distribution of dust 
particles, we applied the following procedure: We computed the dynamical trajectory of 
each particle and kept a record of the particle's orbital elements and positions with a certain 
time interval. These data were then used to characterize the orbital elements and positions 
of many particles over the entire time span, from an initial instant intil the particle's death 
(impact on planet, the Sun, or particle's ejection from the Solar system). For a stationary 
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system, such as the Sun, the constant number of dust grains (being in a balance between 
their production by minor bodies and eventual disappearance in due course of a drift toward 
the Sun under the P-R drag), and no planets included, this approach would not require 
a detailed proof. For a system that incorporates several planets, any particular grain's 
trajectory cannot be considered as a stationary one, the cloning of one dust particle into 
many others as outlined above is justified by the following arguments. 

After a dust particle starts its journey, gravitational scattering on the planets causes 
chaotic change of the particle's orbital elements. For example, we found that the orbital 
parameters of a particle of r = 5-10//m changed by ^ 1% in 10^ yrs and by ^ 10% in 10^ 
yrs. In other words, during the particle's lifetime of 2 ■ 10^ yrs, there are hundreds of strong 
changes in its trajectory. Due just to gravitational scatterings on the giant planets, the 
scattered particles rapidly forget their initial conditions (Levison & Duncan 1997). Unless 
the number of particle trajectories is very small, the computed distribution function of 
dust particles depends only weakly on initial (a, e, i)-orbital elements for the particular 
trajectory, and practically does not depend on the time the trajectory starts. Therefore, 
in deriving the particle distribution function for such a highly chaotic system as the dust 
population in the outer Solar system, the initial positions of kuiperoids as the dust sources 
are much more important than the start time for particles. 

For the present paper, we simulated 280 trajectories of dust particles starting from the 
apocenter and pericenter of 100 KBOs, which produced 1.2 x 10^^ particle positions. In 
the course of our computations of each particle's trajectory, the following procedure was 
applied: 

1. The computed particle positions were sorted on a 3D-grid containing 

45 X 180 X 244 ^2-10^ cells with steps in (heliocentric latitude </?, longitude 
A, and radius R) of (2°, 2°, 0.025 R [AU]). 

2. The computer recorded the particle's (a, e, i) orbital elements and {x, y, z) -coordinates 
once per revolution of Neptune around the Sun (i.e. every 164.8 years). This enabled 
us to create an auxiliary file containing 1.2 • 10^^/6 • 10^ = 2 • 10^ particle positions. 
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These coordinates were sorted into two 2D data files: a 100 x 1000 array in the 
(a, e)-plane (a < 150 AU) and a 180 x 1000 array in the (a, i)-space. The following 
bins were used: Aa = 0.15 AU, Ae = 0.01, Ai = 0.5°. A few auxiliary ID files dealing 
with particle distributions in semimajor axis, n(a), in pericentric distance, n(g), and 
in radius, n(r), all used A = 0.3 AU. 

Let us estimate the number of trajectories which would be sufficient to derive the 
stationary distribution function. Suppose that every cell of the (a, e, i)-space or {x, y, z)- 
space is permeated by N trajectories. The necessary condition to get a robust distribution 
of particles in any cell is that a very large number of trajectories. A'" >> 1, visited that 
cell. For practical reasons, we adopt N ~ 10. Since the trajectories are highly chaotic, it 
does not matter whether a particular cell is visited, say, 10 times by the same particle or 
by 10 different particles just once. We found that a particle of r = 5-10/xm requires about 
2 • 10^ yrs (on average) to change one (a, e)-space cell for another, so that it visits ~ 10^ 
cells during its lifetime. Therefore, on a 2D (a, e)-grid containing 100 x 1000 = 10^ cells, 
one needs to simulate ~ 100 different trajectories to achieve, on average, ~ 10 particle 
visits per cell (actually, this number may be as large as a few hundred in the densest 
regions and ^ 10 in the rarefied ones). For the present paper, we similated 80 trajectories 
of r = 5-10//m particles and 200 trajectories of r = 5-10/xm particles, which provides a 
resonable statictics with a rather low level of noise. 

We employ a second-order imphcit numerical integrator described in Taidakova (1997) 
and Taidakova & Gorkavyi (1999). Its basic features are outlined in the Appendix. As 
shown there, for dust particles produced by typical kuiperoids the dissipative integral of 
motion is conserved with an accuracy of 10~^ — 10~^, which is appropriate to explore such 
a highly chaotic system as dust population in the Solar system. 

5. The Results: Dust Belts and Their Resonant Structure 

Using the above computational approach, we have simulated the spatial structure of 
the IPD cloud between ~ 0.5 and ~ 100 AU, and determined distributions in the phase 
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space of orbital elements. Here we present the results for dust particles of radius 1-2 jj,m 
produced by 100 Kuiper belt objects, and particles of radius 5-10 /im produced by 40 
KBOs. The initial conditions for orbit integration were taken in accordance with equations 
(l)-(2). 

The distributions have been computed with two different values of the P-R parameter 
(the radiation pressure to gravitational force ratio) P = 0.057 and 0.285, and with the solar 
wind drag to the PR-drag ratio =0.35 (Gustafson 1994). Details of computational runs 
are given in Table 1. The total CPU time of those computations with a 450 MHz PC was 
about 2 months. 

5.1. Distribution of Dust in Phase Space 

Distribution of Dust in (a, e)- and (a, i) -space. A representative dynamical trajectory 
of a dust particle on the phase plane is shown in Fig. 1, where the changing role of different 
forces at different parts of the trajectory is clearly seen. The trajectory is ended with 
the ejection of the particle from the Solar system by Jupiter, which is typical for most 
particles. The rest of particles eventually find their death on the the Sun or the planets. 
Our computations indicate that about 11.5% of 1-2/xm particles penetrate the innermost 
zone of the Solar system, whereas 88.5% of particles are ejected. Similarly, 13.8% of 5-10//m 
particles penetrate the Earth zone, whereas 85% are ejected from the Solar system (and 
one of the 80 particles, i.e. 1.2%, fell onto Neptune). This ~ 12% fraction of kuiperoidal 
particles penetrating the inner Solar system is consistent with what was found by Liou, 
Zook, & Dermott (1997). This fraction is not very sensitive to the particle size, which can 
be explained as foUowss: although smaller particles drift more quickly toward the Sun under 
the P-R drag, they are more easily ejected from the Solar system due to the combined 
action of the solar radiative pressure and gravitational scattering by the planets. 

The overall picture of dust distribution in the Solar system obtained by summation of 
the computed trajectories is shown in Figs. 2a,b using 80 trajectories of large (r = 5-10//m) 
particles and in Figs. 2c, d using 200 trajectories of small (r = 1-2/xm) particles. We find 
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the simulated dust distribution be highly non-uniform. All three dynamical classes of dust 
particles listed at the end of Sect. 3 are clearly seen on the (a, e)-plane of Fig. 2c: 

- drift particles with a maximum density at a > 45 AU for l-2//m particles and at a > 30 

AU for 5-10/xm particles; 

- resonant particles producing numerous dense populations in the Neptunian and other 

resonances; the resonant population of large particles has a bigger contrast with the 
background than the population of small particles; 

- scattered particles stretched along the right boundary of the planet crossing zone (i.e. 

the particle pericenters are located near the planet's orbit), especially for Neptune, 
Saturn, and Jupiter. Scattered particles tend to avoid resonant orbits thereby forming 
resonant gaps. The Neptunian region possesses the most dense population of scattered 
particles. 

The dust distribution in the (a, e)-space shown in Figs. 2a,c indicates that each planet 
governs a dust belt that consists of both resonant and scattered particles. If the 

resonant particles dominate in the belt, then the latter is associated with a spatial excess of 
dust particles (like the Neptunian belt described in Sect. 5.2). A different situation arises 
if the belt consists mostly of scattered particles: in this case, the belt is characterized by 
a density minimum (like the Jovian belt). In both cases, the belt is revealed and can be 
seen as a maximum in the particle distribution in distance of pericenter, as will be shown 
at the end of this Section. A criterion for a dust particle to be included into the belt is a 
substantial dynamical interaction of the particle with the host planet, either by resonance 
or by gravitational scattering. 

As can be seen from the dust particle distribution on the (a, i)-plane, the inclinations 
of kuiperoidal particles, on average, substantially increase due to gravitational scattering as 
the particles move from the Kuiper belt toward Jupiter. 

Distribution of Dust Particles in Semimajor Axis and Their Resonant Structure. Fig. 3 
shows the distribution of kuiperoidal dust in semimajor axis. The vertical coordinate is a 
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measure of the number of particles within each 0.3 AU bin per trajectory. Fig. 3 reveals a 
rich resonant structure of each dust belt. Arrows show positions of particular resonances. 
One can see that large-size particles are more easily captured in (or spend more time 
within) the resonances than small-size particles. Furthermore, the smaller the particle size, 
the smaller is the contrast of the resonant structure to the background, which confirms a 
similar conclusion reached by Liou & Zook (1999). Large particles with a > 50 AU are 
mostly scattered particles, therefore one can see resonant gaps in their distribution rather 
than resonant peaks, whereas peaks arc seen at smaller a = 35 — 50 AU. Indeed, as Fig. 3 
demonstrates, the fraction of resonant particles in the Neptunian dust belt is very large. 

The following resonant features seen in Fig. 3 are worth mentioning: 

(i) in the scattered dust component, the gaps at the resonances 5:2, 7:2, 4:1, 5:1, 6:1 etc., 
are pronounced; 

(a) for captured particles, peaks at the resonances 6:5, 4:3, 3:2, 5:3, 7:4, 9:5, 2:1, 4:1, 5:1, 
6:1, etc., are pronounced; 

Like the cometary belts simulated in OGT (2000a,b), the simulated dust belts indicate 
a complex structure containing many families of captured resonant particles and gaps. 
While resonances in the cometary population can form in a non-dissipative way, the resonant 
capture of dust particles occurs dissipatively, and this process takes place both inside and 
outside the crossing zone. For (3 0, resonances are shifted by a factor (1-/3)^/^. The 
larger the value of (3, the larger is the drift velocity and the smaller is the probability of a 
resonant capture. 

In a resonance (j -|- while the eccentricity is close to the maximum, 

Cmax = ^0.4/(j -I- 1) (Weidenschilling & Jackson 1993), the particle's resonant hfetime is 
expected to be long. The resonances seen in Figs. 2a and 2c demonstrate this kind of 
behaviour: the larger is eccentricity of particles in a given resonance, the more abundant 
is their population. This results in a two-hump structure seen, for large-size particle 
distribution, between g f=i 25 and 40 AU in Fig. 4 and between R ^ 25 and 45 AU in Fig. 5. 
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The characteristic shape of such dust structures follows from the Kepler laws (Kessler 1981, 
GOMT 1997b): The inner and outer edges of each structure are given by Uresi^ ~ ^max.) and 
Qres(l + Cmax), respectively. For all resonances, in a good approximation, the position of the 
inner edge is ~ 0.85apianet- At this position, we expect to find rather sharp inner edges and 
steps in the dust density distribution (somewhat depending on the particle size) for each 
giant planet's dust belt. 

Direct measurements of the expected dust density maximum at ~ 27 — 29 AU by 
dust detectors on spacecraft would be a strong confirmation of the predicted Neptunian 
dust belt, as well as its resonant nature simulated in the present work. 

Distribution of Dust Particles in Pericenter Distance. Fig. 4 demonstrates that the 
pericentric distances of many kuiperoidal particles are located close to the orbit of each 
giant planet (see also Figs. 2a,c). One important difference between the resonant and 
the scattered components of a belt is in the value of distance of pericenter, which is 
Qr ~ (0.85 — 1.0)api for the resonant component (see the location of the Neptunian belt 
relative to the Neptune's orbit shown as N in Fig. 4), and ^ (1.0 — l.l)ap; for the 
scattered component (see the locations of the Jovian and Saturnian belts relative to the 
respective planets shown as J and S). Thus the distribution in pericentric distance provides 
an additional strong argument in favor of the four dust belts. Earlier (OGT 2000a,b), 
we have found that cometary bodies concentrate into belts near each giant planet's orbit. 
As the parameter /? decreases with increase of particle size, the efficiency of gravitational 
scattering increases, which makes the pericentric distribution of large dust grains like that 
of comets. 

5.2. Spatial Distribution of Kuiperoidal Dust Particles 

Fig. 5 presents the column density of dust particles [using all recorded particle positions 
{x, y, z)] as a function of heliocentric distance. Going inward, the following features of dust 
density distribution are worth emphasizing: 
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(i) an approximately constant column density of dust (as well as volume density) between 

50 and 10 AU; 

(ii) a sharp decrease of both column and volume density of dust between 10 and 4 AU, 

which is due to the ejection of particles by Jupiter and Saturn; 

(in) an increase of column density of dust (accompanying by an even steeper increase of 
volume density near the ecliptic plane) at it! < 4 AU. 

The column density of kuiperoidal dust forms a plateau between 50 and 10 AU, and 
this seems to be the most remarkable result. For 5-10 /xm dust, the value of column density 
is several times higher than that for 1-2 yum grains, which is explained by a slower rate of 
evolution of larger size particles, mostly due to the P-R drag. The influence of gravitational 
scattering and resonances is also different for particles of different sizes, so that the shape 
of the plateau, as can be seen from Fig. 5, depends upon particle size. 

The distribution of the number density of dust in the ecliptic plane, which is of obvious 
practical interest to interpret the data of space missions, is qualitatively similar to column 
density graph. This is because the particle inclination, on average, decreases outward (see 
Figs. 2b, d). A region of an elevated density of kuiperoidal dust, which is associated with 
the Neptune's orbits, can be seen in Fig. 5. In the ecliptic plane, the Neptunian dust belt is 
expected to have the largest number density of particles. The major part of the simulated 
Neptune's dust belt consisting mostly of resonant particles is located between 25 and 45 
AU and forms a flat dense disk. The simulated Uranian, Saturnian, and Jovian dust belts 
basically overlap and form complex structures which, in their central parts, are less dense 
than the Neptunian dust belt. 

There is a remarkable density minimum between Mars and Jupiter. This minimum 
is due to the fact that Jupiter either ejects from the Solar system or transfers to more 
inclined and eccentric orbits an appreciable part of the dust drifting toward the Sun. An 
increase of dust number density going inward from Mars to Earth can be explained by the 
dominant role of the P-R drag here, which would result in the number density distribution 
n{R) oc for circular orbits or n{R) oc i?"*^^"^^ for eccentric orbits (see GOMT 1997b 
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and refs. therein). 

Fig. 6 shows a 2D section of the spatial structure of kuiperoidal dust number density 
in the region up to 10 AU perpendicular to the ecliptic plane (edge-on view). An important 
feature found in our simulations is the existence of a new quasi-stationary, highly inclined 
dust population with pericenters near Jupiter and Saturn. This population is seen in 
Fig. 6 as a sharp increase in dust density beyond the Jupiter orbit, which looks like 
a 'Chinese wall'. This structure is found to be steeper and denser for larger particles. 
The quasi-stationarity of the structure results from a balance between the tendencies for 
particle's semimajor axis a and eccentricity e to increase due to gravitational scattering 
on the planet and to decrease due to the P-R drag. The particle inclinations increase 
substantially due to gravitational scattering and resonances. 



6. Discussion and Conclusions 

We have explored whether the 280 particle trajectories used in the present study 
are sufficient to provide reliable results. As can be seen from Fig. 5, the use of just 20 
particle trajectories already reveals the major features of large-size particle distribution, 
but for increased accuracy we used 80 trajectories. To explore the detailed distribution 
of kuiperoidal dust in phase space, we need as many as 80 trajectories to reveal the basic 
resonant features (see Figs. 2a,b). The comparison of Figs. 2a,b with Figs. 2c, d (where 
200 trajectories are used) indicates that a larger number of computed trajectories leads 
to a decrease in random fluctuations; gives a clearer picture of strong resonances free of 
discontinuities; reveals a larger number of weak resonances; and results in a more reliable 
picture of resonant gaps and near-resonant accumulations in the scattered component of 
dust. Although even larger statistics would certainly provide further improvements, we do 
not think that it would qualitatively change the results. 

This study employs a number of approximations and simplifications (circular planetary 
orbits with zero inclinations, neglecting other components of dust, including the interstellar 
dust), which will be treated more accurately in further work. Flynn (1996) and Liou et al. 
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(1997) claim that large (^9 fim.) kuiperoidal dust is destroyed by collisions with interstellar 
dust. Nevertheless, it is instructive to compare our results with available observational 
data. 

Pioneers 10 and 11 as well as Voyagers 1 and 2 detected a large number of dust 
particles in the region between Jupiter and Neptune (Humes 1980; Gurnett et al. 1997). 
These data show that the number density of dust particles in the outer Solar system is 
approximately distance- independent. As can be seen in Fig. 5, the simulated distribution 
of kuiperoidal dust can explain this constancy. This qualitative result will be quantified in 
our further work, after improving our modelling and taking into account relative velocities 
of dust particles and the spacecraft. 

The approximate constancy of the dust density between 10 and 50 AU found in the our 
simulations differentiates the kuiperoidal dust from both the asteroidal and cometary dust, 
whose number densities are known to appreciably decrease with distance from the Sun (for 
interpretation, see GOMT 1997b and refs. therein). 

Our simulations offer a high-resolution, 3-D model of the kuiperoidal cloud on a grid 
of 2 million cells containing 115 billion computed positions of dust particles. The major 
conclusion reached in the present simulations are as follows: 

1. The simulated dust distribution is highly non-uniform. Moving inward, the column 
and volume density of kuiperoidal dust is approximately constant at heliocentric distances 
from 50 to 10 AU, sharply decreases between 10 and 4 AU giving a deep minimum, after 
which the column density of dust increases (accompanied by a constant or decreasing 
volume density near the ecliptic plane) at i? < 4 AU. 

2. We find a new quasi-stationary, highly inclined dust population with pericenters near 
all giant planets. This quasi-stationarity results from a balance between the tendencies for 
the particle semimajor axis a and eccentricity e to increase due to gravitational scattering 
on the planet and to decrease due to the P-R drag. The particle inclinations i increase 
substantially due to gravitational perturbations from Jupiter and Saturn. 

3. Most of the dust is concentrated into four belts consisting of resonant and scattered 
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particles associated with the orbits of the four giant planets. These belts are chiefly of 
dynamical nature. Some of them, such as the Neptunian belt, disclose first of all as dust 
structures. The others have very distinct distributions in the phase space: for instance, the 
Jovian (Saturnian) belt is characterized by a substantial excess of highly eccentric particles 
with the maximum of dust distribution in pericentric distance at q a; 5 (10) AU. The 
Saturnian belt could be discriminated by an excess of dust particles in the resonances 4:3, 
3:2, and 2:1. 

4. The simulated dust belts reveal a complex resonant structure containing many 
families of gaps and resonant maxima. The particles are either dissipatively captured into 
exteriour resonances (usually outside the crossing zone) or form gaps (usually inside the 
crossing zone). 

5. A rather long life time in each resonance, while the eccentricity is close to the 
maximal one, results in a steep rise of dust density at the innermost edge of the resonant 
component in all dust belts, especially the Neptunian one. Rather sharp inner edges and 
'steps' in the dust density distribution are expected to characterize each giant planet's dust 
belt at 0.85api, where Opi is the planet's semimajor axis. Neptune's dust belt is expected to 
have both the largest 'step' and number density of particles in the echptic plane (Fig. 5). 
Direct detection of a dust density maximum at fti 27 — 29 AU in Neptune's zone would 
test the simulated dust distribution. The Neptunian and the other belts would be a 
challenging target to discover and measure by space missions. 

The resonant features of dust distributions near giant planets can serve as signatures 
of exo- planets in the circumstellar disks (Ozernoy et al. 2000c, Gorkavyi et al. 2000a). The 
kuiperoidal dust is likely to be a contributor of the zodiacal light emission in the Solar 
system, which is analyzed in Gorkavyi et al. (2000b). 
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Appendix 



To simulate tlie dissipationlcss dynamics of comets as well as dissipativc dynamics of 
dust particles orbiting around a star, we use an implicit second-order integrator (Potter, 
1973; Taidakova, 1997; Taidakova & Gorkavyi, 1999; Pridman & Gorkavyi, 1999) in a 
rotating (comoving with Neptune) coordinate system. The latter is convenient to show the 
planetary resonant structure as a stationary pattern. 

The equations of motion of a dust particle in the gravitational field of the Sun and the 
planets in this coordinate system take the form (Taidakova 1990, 1997): 



X — 2y -\- X -\- Fj. 
y = -2x + y + Fy 



where F^, Fy, F^ are the components of the sum of the gravitational forces and the 
PR-drag. Integration employs the following equations (Taidakova 1990, 1997): 
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where x'^+^l = + v^^^At/2 ; yl'^+^l = + v^^^At/2 ; ^I'^+^l = + vMAi/2 ; and 
Fi:t^l = F ^["+^1, ^["+^1) . 



We have tested our integrator for a non-conservative system that includes the Sun 
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and a test particle and is governed by the Poynting-Robertson drag (two values of the 
parameter (3 — 0.285 and 0.057 have been used). As initial conditions for {a, e, i}, we used 
the following sets: {39 AU, 0.25, 10°} and {45 AU, 0, 10°}. The integration time step 
was taken in the range 16 to 160 days (it was smaller, the closer the test particle to the 
star). Fig. 7 shows the accuracy of our integrator evaluated by the change of the first 
dissipative integral of motion with time, 5C — C{t) — C(0), where C — ^ (see, e.g., 
GOMT 1997b). It can be seen that, as the particle approaches the star, the integration 
error increases, but it never exceeds 1% and is much better during the larger part of the 
trajectory. We emphasize that, since in the real Solar system any dust particle experiences 
strong gravitational perturbations from the planets so that the particle's trajectory is highly 
chaotic, an accumulation of dissipative integral errors as small as shown in Fig. 7 is of no 
importance. 

Thus, our integrator takes into account close approaches with planets, which occur 
frequently for dust particles spiraling toward the star due to the PR-drag, and demonstrates 
good results in terms of stability of the integration error. In addition, this integrator is 
1.5 — 1.9 times faster than the ordinary 2nd-order Runge-Kutta integrator (see Taidakova 
1997). Although our integrator is not as speedy as symplectic integrators, it has an 
important advantage since it is applicable to dissipative systems, along with non-dissipative 
ones. 

We found that our integrator employed in a non-rotating reference frame has a similar 
accuracy but is several times faster than in a rotating system (Taidakova et al. 2000, in 
preparation). As an example of the use of our integrator in a non-rotating frame, we have 
recently simulated a warp observed in the circumstellar disk of Beta Pictoris (Gorkavyi et 
al. 2000c). 
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Details of Computational Runs 



Particle lifetime Number of recorded Number of computed 

(in Myrs) (^) (a, e, i)-elements(^) spatial positions 

(in 10^) (in 10^°) 

Grains of r = 5-10/xm 
No planets, 

20 grains used 10.1 1.2 0.7 
Planets included, 

80 grains used 20.3 9.9 5.5 

Grains of r — l-2//m 
No planets, 

20 grains used 4.1 0.5 0.3 
Planets included, 

200 grains used 6.7 8.1 6.0 



(^) until the particle impacts the Sun, or a planet, or is ejected from the Solar system. An 
average value of the lifetime is given. A larger value of the lifetime for larger-size 
grains is explained by a slower PR-drag. Presence of planets additionally increases 
the particle lifetime due to captures into resonances. 

(^) taken with the time step = 1 Neptune's revolution about the Sun. 
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Figure Captions 

Figure 1. Representative trajectory of a dust particle of of r = 5-10/xm (/3 = 0.057) 
on the planes of orbital coordinates {a, e} (panel a) and {a, i} (panel b). The trajectory 
presents the particle positions taken every 5 x 10^ yrs. Dashed curves show the boundaries 
of the crossing zones of the four giant planets. Diamond indicates the particle's initial 
position. Numbers i to 5 mark the dominating dynamical process on the given part of 
the trajectory: 1 - a, 'jump' from the parent body (KBO) due to the solar pressure (see 
Eqs. 1 and 2); 2 - drift of particles due to the P-R drag; 3 - resonant capture into the 3:2 
resonance with Neptune, which results in a balance between the P-R drag and gravitational 
influence of the planet; ^-5 - gravitational scattering of the particle by giant planets, with 
eventual ejection out of the system by Jupiter. 

Figure 2. 2D density of the kuiperoidal dust on the plane of orbital coordinates, 
with bins Aa = 0.15 AU, Ae = 0.01, Ai = 0.5°. To represent the number of particles in 
each cell, a decimal-logarithm grey scale is employed, i.e. each shade differs 10-fold from 
the neighboring one. Numerous resonant lines and gaps are seen. The boundaries of the 
crossing zones of the four giant planets are indicated by solid curves. Positions of the 
first 40 (a, b) or 100 (c,d) Kuiper belt objects taken from Marsden (1998) are shown by 
diamonds. Numerous resonant structures are clearly seen. 

a, b. 80 dust particles of r = 5-10/xm (/3 = 0.057) start their journey from the 
pericenters and the apocenters of orbits of 40 kuiperoids. 

c, d. 200 dust particles of radius r = 1-2/xm (/3 = 0.285) start their journey from the 
pericenters and the apocenters of orbits of 100 kuiperoids. 

Figure 3. Distribution of kuiperoidal dust in semimajor axis, in terms of the number of 
particles per bin, n{a), averaged over each trajectory. The bin size Aa = 0.01 aNeptune = 0.3 
AU. The distributions of small (r = 1-2/xm) and large (r = 5-10/im) dust particles are 
shown by dashed and solid lines, respectively. Various resonant structures are indicated by 
arrows, which are heavy for large particles and thin for small ones (J, S, and U stand for 
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the Jovian, Saturnian, and Uranian resonances, respectively; all other resonances are with 
Neptune) . 

Figure 4. Distribution of kuiperoidal dust in the distance of pericenter, in terms of 
the number of particles per bin, n{a), averaged over each trajectory. The bin size Aa = 0.3 
AU. The distributions of small (r = 1-2/xm) and large (r = 5-10/im) dust particles are 
shown by dashed and solid lines, respectively. 

Figure 5. Column density of kuiperoidal dust population, in terms of number of 
particles per 0.3 AU bin averaged over one trajectory, as a function of heliocentric distance 
(solid line). The distributions of small (r = 1-2/im) and large (r = 5-10/im) dust particles 
are shown by thin and heavy lines, respectively. Dotted line shows a distribution of large 
particles obtained with the use of only 25% (i.e. 20) of the available particle trajectories, 
which indicates that the results depend rather weakly upon the number of trajectories used. 
Dashed and dashed-dottcd lines indicate the surface density of kuiperoidal dust computed 
from 20 particle trajectories in the absence of planets. Larger-size particles, due to a slower 
motion, form a denser dust population. 

Figure 6. Density profile of a small, r — l-2//m, kuiperoidal dust perpendicular to the 
ecliptic plane, N[R, Z), at heliocentric distances up to 10 AU. To represent the number of 
particles in each cell, a natural-logarithm grey scale is used, i.e. each shade differs e-fold 
from the neighboring one. A remarkable density minimum between Mars and Jupiter is 
clearly seen. 

Figure 7. An average numerical error in the dissipative integral of motion along the 
particle trajectory (shown as a function of semimajor axis) for a dust particle drifting 
under the P-R drag toward the Sun. Trajectories starting from the resonant kuiperoids are 
marked with 1 and 3, while those starting from the fiat component of KBOs are marked 
with 2 and 4- Initial orbital positions {gq, cq, io} of the parent KBOs are as follows: 

la,lp - {39 AU, 0.25, 10°}, start from apocenter and pericenter, respectively; /3=0.285. 
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2 - {45 AU, 0, i = 10°}; /3=0.285. 

3a,3p - {39 AU, 0.25, i — 10°}, start from apocenter and pericenter, respectively; 
/5=0.057. 

4 - {45 AU, 0, i = 10°}; /3=0.057. 

For an overwhelming majority of trajectories, the integration errors, 10~^ — 10~^, are 
well within the acceptable limits. 
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